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One possibility for explaining the apparent accelerating expansion of the universe is that we 
live in the center of a spherically inhomogeneous universe. Although current observations cannot 
fully distinguish ACDM and these inhomogeneous models, direct measurement of the acceleration 
of the universe can be a powerful tool in probing them. We have shown that, if ACDM is the 
correct model, DECIGO/BBO would be able to detect the positive redshift drift (which is the 
time evolution of the source redshift z) in 3-5 year gravitational wave (GW) observations from 
neutron-star binaries, which enables us to rule out any Lemaitre-Tolman-Bondi (LTB) void model 
with monotonically increasing density profile. We may even be able to rule out any LTB model 
unless we allow unrealistically steep density profile at z ~ 0. This test can be performed with GW 
observations alone, without any reference to electromagnetic observations, and is more powerful 
than the redshift drift measurement using Lyman a forest. 
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Introduction: If we assume that our universe is homo- 
geneous and isotropic, current cosmological observations 
(e.g. type la supernovae (SNe) indicate the acceler- 
ating expansion of the universe. Dark energy and mod- 
ification of the gravitational theory are two candidates 
that can explain these observational results. However, 
once we allow for the possibility that our universe has 
cosmological-scale spherical inhomogeneity, with the ob- 
server at the center, the observations can be explained 
without introducing the unknown dark energy or alter- 
native theories of gravity. In such models, the Coperni- 
can Principle is apparently violated and the expansion of 
the universe is not necessarily accelerating. As such, the 
direct detection of the cosmic acceleration is very useful 
in testing them. Although there are several proposals 
(e.g. Refs. @, 01) for the direct detection of the cosmic 
acceleration, it has not been observed yet. The direct 
detection of the acceleration of the universe provides not 
only a key to solving the dark energy problem, which 
is one of the biggest challenges in cosmology, but also a 
critical test of the Copernican Principle, which is one of 
the most essential "principles" in cosmology. 

The simplest example of the inhomogeneous model is 
the Lemaitre-Tolman-Bondi (LTB) spacetime, which is 
a spherically symmetric, dust solution of the Einstein 
Equations. If we live at the center of the LTB spacetime 
with a Gpc-scale void, the apparent acceleration of the 
universe can be explained. The LTB model has been par- 
tially tested with cosmological observations such as, the 
cosmic microwave background, baryon acoustic oscilla- 
tions, the kinetic Sunyacv-Zcldovich effect etc., but it has 
not been completely ruled out yet (see e.g. Ref. [J). One 
of the fundamental difficulties in testing the LTB mod- 
els is due to the ambiguities in the primordial spectrum. 
In the standard cosmology, we mostly assume a scale 
invariant spectrum based on the inflationary paradigm 
that makes our universe homogeneous in the early epoch. 



However, once we introduce cosmological scale inhomo- 
geneities and violate the Copernican Principle, there is no 
strong motivation to assume inflation in the early epoch. 
Therefore, observations which are not affected by the pri- 
mordial spectrum are very crucial in testing a wide class 
of LTB models. One such observation is the redshift- 
distance relation of type-la SNe. However, it is known 
that one can construct a LTB void model that exactly 
reproduces the redshift-distance relation of the ACDM 
model m. Therefore, we need other observations that do 
not depend on the primordial information. 

Redshift-drift measurement is the one that meets our 
demands @ (see also Refs. 0, B)- Redshift drift is 
the time evolution of the redshift due to the accelera- 
tion of the universe, hence its detection means the direct 
measurement of the cosmic acceleration. 

In the Friedmann-Lemaitre-Robcrtson- Walker 
(FLRW) spacetime, the redshift drift is given as 

AfZ ~ HoAto (^l + z — I where Ato denotes the 

observation period, and Hq and II{z) are the Hubble 
parameter at present and at redshift z 0, respectively. 
In the ACDM universe, AfZ is positive in the range 
z = 0-2 [3] . On the other hand, typical LTB models that 
can explain observations usually have negative Atz at 
any z (see e.g. Ref. Q)- Recently, Yoo et al. [1| showed 
that if the matter density monotonically increases with 
the radial distance r, A-tZ must be negative for any 
z > 0. Furthermore, they have shown that for any 
LTB density profile, dAtz/dz < for z < 1. This 
indicates that, unlike the redshift-distance relation, it 
is impossible to exactly reproduce AtZ of the ACDM 
model with LTB models. Therefore, it is crucial to 
measure the sign of Atz at z < 2 (especially z <C 1) m 
order to distinguish the ACDM and LTB models. 

The order of magnitude of AtZ is roughly given as 
the cosmic age divided by the observation time, hence 
AtZ ^ 10~^° X (obs. time)/lyr. It is this tiny value that 
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makes it difficult to measure Ajz witfi current tectmol- 
ogy. Recently, Quartin and Amendola have shown 
that by measuring the shift of the Lyman a forest of 
quasar spectrum at z = 2~5 with the proposed E-ELT 
instrument CODEX for 10 yrs, it will be possible to 
distinguish ACDM and typical LTB void models. How- 
ever, CODEX would not be able to measure Atz at low z 
since the Lyman a forest can be measured from ground 
only at z > 1.7 They can only test typical LTB 
models and not generic ones. 

In this letter, we estimate how accurately we can mea- 
sure the redshift drift with future gravitational wave 
(GW) interferometers. It seems that DECIGO d [l^ 
and BBO [ll| are the only proposed detectors that can 
measure AfZ at z < 2. These are complementary obser- 
vations to electromagnetic ones. We consider neutron- 
star (NS) binaries as GW sources, which are often called 
standard sirens and can be unique tools to probe the 
cosmic expansion [l^ - [l^ . When the expansion is accel- 
erating, we may find an additional phase shift in gravi- 
tational waveforms associated with the redshift drift 
We assume that ACDM is the correct model and estimate 
whether we can tell the positivity of the redshift drift at 
low z with GW observations. 

Throughout this letter, we use the unit G = c = 1. 

Measuring the redshift drift with GWs: Let us first de- 
rive the correction to the GW phase due to the accelerat- 
ing expansion of the universe. We here consider a binary 
consisting of two bodies with masses mi and TO2 . We de- 
fine h{At) as the observed GW where At = tc — t denotes 
the time to coalescence measured in the observer frame 
with tc representing the coalescence time. The Fourier 
component of this waveform is written as 

/oo 
dt e^'''f%{At) 
-OO 

/OO 
dAt e-2"/'^*/i(A<) . (1) 
-OO 

Now, At can be related to AT = (1 + zJAie as 

At = AT + A(zc)Ar^ (2) 

where Zc is the source redshift at coalescence. Ate is the 
time to coalescence measured in the source frame and 
X{z) is the acceleration parameter defined as X{z) = 

^{l- (1+1)^0 )■ Notice that X{z) is related to the 
redshift drift Atz as Ajz = 2(1 ^ z)AtoX{z) in FLRW 
spacetime. Then, h{At) can be re-expressed as a func- 
tion of AT(At) as h{At) = H{AT{At)), where H{AT) 
corresponds to GWs without cosmic acceleration. By 
substituting this equation and Eq. ([2|) into Eq. ([T]), we 
get 

/OO 
dAT' e-^'^'f^^' H{AT') 
'OO 



By using the stationary phase approximation [16|, 
g-2irj/x(zc)AT integrand can be pulled out 

of the integral with AT' replaced by Ar(/) = 
b{%TTM:,f)-^I^M^ 0, where = Mil + Zc)rf'^ de- 
notes the redshifted chirp mass with M = (mi -I- m2) 
and rj = mim2/M^ representing the total mass and the 
symmetric mass ratio, respectively. Then, we obtain 
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= »*aoc(/) 



h{I)\ 



no acccl 



where 



*acc(/) 



2^/X(z,)Ar(/)2 

{f)-^X{z,)M.x-\, 



768 



(4) 



(5) 



with X = {-KM.Sf/'^ and *Ar(/) = ^(7rX^/)-5/3^ and 



no acccl 



/OO 
dAT' e-^'^'i'^'^' H{AT') , 
-OO 



(6) 

which corresponds to the gravitational waveform in the 
Fourier domain without cosmic acceleration. The above 
^'acc(/) matches with the ones found in Refs. 0, [HI. A 
term with a;" represents the rt-th post-Newtonian (FN) 
order relative to the leading Newtonian phase \I'Ar(/), 
hence this is a "-4PN" correction. For h{f)\ ,, 
we use the restricted- 2PN waveform including spin- 
orbit coupling at 1.5PN order [3l| which is given in 
Eq. (26) jsil of Ref. [13, where "restricted" means that 
we only take the leading Newtonian quadrupole contri- 
bution to the amplitude and neglect contributions from 
higher harmonics. 

Numerical Setups: The measurement accuracy on 
Atz has been estimated mainly using ultimate DECIGO 
(which is three orders of magnitude more sensitive than 
DECIGO) I, Il5l |. In this letter, we improve their analy- 
ses in the following way and apply the result to probe 
the inhomogcneity of the universe: (I) We use DE- 
CIGO/BBO, including the confusion noises from white 
dwarf (WD) binaries. (II) Rather than sky-averaged 
analysis, we perform Monte Carlo simulations by ran- 
domly distributing the directions and orientations of 
sources. (Ill) We include the spin-orbit coupling into 
binary parameters. (IV) We use the merger rate that re- 
flects the star formation history rather than the redshift 
independent rate. (V) We do not require that the source 
redshift z is known by identifying the host galaxies with 
EM observations. (VI) Since we arc only interested in the 
positivity of the redshift drift, we consider the cumulative 
acceleration parameter x'^^^'^\z) defined as 



Ar'"^(z,)^5:AH(z„), 



(7) 



xe 



-27ri/X(ze)AT'^ 



(3) 



where Xh{z) = X{z)/Hq and z„ = (n — 0.5)6z with Sz 
representing the size of the redshift bin. 
Here, we take the binary parameters as 

9' - {lnM,,lnr],p,tc,(l)c,DL,9s,^s,0L,^L,XH) ■ (8) 
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/9 is the spin-orbit couphng parameter while <f>c repre- 
sents the coalescence phase and Di^ denotes the luminos- 
ity distance. (^s,<^s) give the direction of the source in 
the barycentric frame which is tied to the echptic and 
centered in the solar system barycenter (see Fig. 1 of 
Ref. and (^l,^l) specify the orientation of the 

source orbital axis in the barycentric frame. 

We estimate the measurement accuracies of binary pa- 
rameters 0^ using Fisher analysis. Assuming that the 
detector noise is stationary and Gaussian, the measure- 
ment accuracy is given as = ^ (r^-^)jj [13 • Here, the 
Fisher matrix F.y is defined as F,,, = {dih\djh), where the 
inner product is defined as 



{A\B) = 4Re 



/tin 



df 



A*U)B{f) 

SnU) ' 



(9) 



with Sn{f) denoting the noise spectrum of a single in- 
terferometer and /in and /fin representing the initial and 
final frequencies of observation. The measurement error 



of xj^"™''(zfc) can be estimated as 



-I 1/2 



(10) 



We use the instrumental noise spectrum of BBO |33j 
and also consider WD/ WD confusion noises which de- 
teriorate the sensitivity below / ^ 0.2Hz. In this let- 
ter, we are interested in the phase shift due to the red- 
shift drift, which appears at relative "-4PN" order. The 
negative PN order indicates that this effect is larger 
when GW frequency is lower. Since the WD/ WD confu- 
sion noise would mask the lower frequency part of DE- 
CIGO/BBO observation window, this effect cannot be 
neglected. We choose /in and /fin in Eq. ([9]) as /i„ = 
(256/5)-3/87r-i7WJ^/*Ato^^^ and /fi„ lOOHz, which 
correspond to the frequency at Ato before coalescence 
and the higher cutoff frequency of the detector, respec- 
tively [13. For fiducial values, we set mi = 7712 = IAMq 
and take tc = 4>c = [3 = 0- We assume that a flat 
ACDM is the correct model and set the fiducial values 
as Hq = 70km/s/Mpc and the cosmological parameters 
as ilrii ~ 0.3 and flA — 0.7. 

The number of binaries AN{z) that exists in each 
redshift bin with the size Sz is estimated as AN{z) — 
Att [aor{z)]^ n{z){dT/dz)SzAto [T9| . where aor{z) ~ 

dz'/H{z') and dr/dz = {{1 + z)H{z)}-^ with H{z) ee 

Hay^Qmi^ + z)^ + Oa. ao, r(z), and r represent cur- 
rent scale factor, comoving distance and proper look back 
time, respectively. ri{z) = riQR{z) shows the NS/NS 
merger rate per unit comoving volume per unit proper 
time, where we assume the merger rate today as no = 



10-6 Mpc~3 yr-i l20| and 



R{z) 




(^<1) 
(1 < z < 5) 
{z>5). 



(11) 



The merger rate evolution against z reflects the current 
observation of star formation history 21 1. 

For each flducial redshift Zn with the redshift bin size 
set as Sz = 0.1, following Refs. [13, HMIl^i '^^ randomly 
generate 10^ sets of (^s, ^s, ^l, 0l) and for each set, we 
calculate ^(F-^).-. Then, we take the average to yield 
[a/ (r~"'")iij jj^jj- The measurement accuracy at each 2„ is 
estimated as 



A9^ = Nrl'^ANizr, 



^1/2 



(F-) 



1/2- 



(12) 



where A^int = 8 represents the effective number of in- 
terferometers [H. DECIGO/BBO has four clusters of 
triangular detectors (sec e.g. Fig. 2 of Ref. [l^ for the 
proposed configurations of DECIGO/BBO) and for sim- 
plicity, we assume that all the clusters are placed on the 
same site. First, we show a rough estimate of the mea- 
surement accuracy of Xh using DECIGO/BBO and then 
show our numerical results. 

Rough Estimate: Let us assume that we observe 
AN{z = 1) = 2.4 X 10^ NS/NS binaries for 5yr obser- 
vations. The measurement accuracy of Xh{z) is mostly 
determined at the lowest frequency of the observation, 
i.e. at / = /syr = 0.073Hz since the effect due to 
the cosmic acceleration is larger on GW signals with 
lower frequency. We first define the squared signal-to- 
noise ratio (SNR) as SNR^ = jj'^-- p^{f)d\nf, where 

p(/)^ = 4A^int/|/i(/)P/5'„(/) expresses the contribution 
to SNR^ at each frequency /. Typical values for p{fin) 
and Xh at this redshift would be p{fin) ~ 1-6 and 
Xh{z = 1) « 0.05, respectively. Roughly speaking, if 
*acc(/)p(/in)VAA^ exceeds 0(1), the effect due to the 
cosmic acceleration can be detected (when we assume 
that there is no correlation between AtZ and other pa- 
rameters). Therefore, we have the measurement accuracy 
of AXh(z = 1) « 0.05 which is comparable to the fidu- 
cial value of Xh{z = 1). Hence, it is marginal whether 
the positivity of the redshift drift can be measured or 
not. We expect that the measurement accuracy would 
be improved when we consider X^^™\zk). 

Numerical Results: In Fig.[l] the (blue) thin line shows 
Xh and its measurement accuracy for 5yr observations 
with DECIGO/BBO. In GW observations, since the red- 
shift z is degenerate with mass parameters [13] , we show 
the result against D^, with corresponding z just for ref- 
erence. We emphasize that the test can be performed 
with GW observations alone by using the Dl~Xh re- 
lation. It can be seen that at z 0.5, DECIGO and 
BBO can marginally probe the positivity of the redshift 
drift. Our result is consistent with the rough estimate 
explained in the previous paragraph. This test can be 
improved, shown by the (red) thick lines, by considering 

the cumulative acceleration parameter a]^"™'. By using 
this quantity, we can probe the positivity of the redshift 
drift with sufficient accuracy. 

Figure [3] shows the ratio between the fiducial value of 
^jcum) measurement error. We can see that 
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FIG. 1; The acceleration parameter Xh (blue thin), the cu- 
mulative acceleration parameter (red thick) and their 
measurement accuracies under ACDM model with 5yr ob- 
servations using DECIGO/BBO. On the horizontal axes, we 
show both Dl and z. The former can be measured from GW 
observations while the latter is degenerate with mass param- 
eters. 
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FIG. 2; The ratio between the fiducial values of X^^"™' and 
its measurement accuracies for lyr (green dotted), Syr (blue 
dashed) and Syr (red solid) observations. This shows to what 
cr-level we can detect the positivity of the redshift drift with. 



it is possible to measure the positivity of Ajz with l.S-cr 
for 3yr observations and 4-cr for Syr observations. Fur- 
thermore, with Syr observations, it may even be possible 
to probe the positivity of Ajz at around z w 0.2 with 2-cr 
confidence. 

Conclusions and Discussions: In this letter, we have 
estimated how accurately we can directly measure the ac- 
celeration of the universe with future space-borne GW in- 
terferometers such as DECIGO and BBO. Assuming that 
the ACDM model is correct, we have shown that, with 3- 
Syr observations, we will be able to measure the positiv- 
ity of the redshift drift, which enables us to rule out any 
LTB void model with a monotonically increasing den- 



sity profile. Also, with Syr observations, DECIGO/BBO 
can measure the positivity of the redshift drift at even 
z ^ 0.2. This indicates that future GW observations will 
be able to rule out any LTB void model unless we allow 
unrealistically steep density gradient at z ~0. Our results 
regarding the use of DECIGO/BBO are complementary 
to the ones proposed by Quartin and Amendola where 
EM observations are used in the sense that these observa- 
tions target different sources at different redshift ranges. 
However, GW observations can be more powerful since 
they can offer a more generic test of the LTB model. 

A variety of dark energy models predict similar values 
of Atz to the ACDM model (see Quartin and Amen- 
dola and references therein). Therefore, we expect 
that the test can be performed with almost the same ac- 
curacy even if the fiducial model is another dark energy 
model. 

Unfortunately, it seems very difficult to measure AtZ 
with ground-based detectors, even with third-generation 
ones like Einstein Telescope (ET) The advantages 

of using DECIGO/BBO over ground-based ones are (i) 
they have a larger number of GW cycles, (ii) the effect 
of AtZ is "-4PN", meaning that we can perform our test 
better with lower- frequency GWs, (iii) they have a longer 
observation time for each binary, and (iv) they detect a 
larger number of NS binaries. 

We need to comment about the peculiar acceleration 
of each binary source. This acts as an additional "noise" 
when measuring the redshift drift. Amendola et al. (25j 
have estimated the peculiar accelerations for typical clus- 
ters and galaxies, and found that they are almost the 
same magnitude as the cosmological acceleration. How- 
ever, they are much smaller than the measurement errors 
of AtZ from a single binary source due to detector noises, 
hence we can safely neglect the effect of peculiar acceler- 
ations (see also Uzan et al. [1^). 

In this letter, we have assumed that the binary orbits 
are circular. If we include the eccentricity, this may be 
degenerate with Ajz since both have large effects when 
the binary separation is large. Thus, in future work, we 
need to estimate how accurately we can determine the 
positivity of AtZ in eccentric binaries. Also, we have 
used a]^"™^ as our estimator because of its simple form, 
but it may not be the best one. It would be interesting 
to find the corresponding one that performs better in 
probing the positivity of the redshift drift. Furthermore, 
we note here that we have used the LTB metric as a 
simple effective model. To be more realistic, we need to 
use more sophisticated models such as the Swiss-Cheese 
model H^l (see a recent discussion by Celerier [H]). We 
leave these issues for future work. 
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